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Abstract 

The cross section for radiative capture of neutron on carbon-14 is calculated using the model- 
independent formalism of halo effective field theory. The dominant contribution from El transition 
is considered, and the cross section is expressed in terms of elastic scattering parameters of the 
effective range expansion. Contributions from both resonant and non-resonant interaction are 
calculated. Significant interference between these leads to a capture contribution that deviates 
from simple Breit-Wigner resonance form. 
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I. INTRODUCTION 



The radiative capture of neutron on carbon-14 14 C(n, 7) 15 C plays an important role in 
astrophysics. It is part of the neutron induced carbon-nitrogen-oxygen (CNO) cycle in the 
helium burning layer of asymptotic giant branch stars and in the core helium burning of 
massive stars [TJ. These neutron induced reactions can lead to appreciable changes in the 
CNO abundances. 14 C(n, 7) 15 C is the slowest reaction in the cycle and leads to substantial 
enrichment of 14 C abundance [1]. In astrophysical scenarios involving inhomogeneous Big 
Bang Nucleosynthesis, the slow 14 C(n, 7) 15 C reaction acts as a bottle neck in the production 
of heavier nuclei A > 14 [21 E]- The 14 C(n, 7) 15 C cross section has been measured in 
direct capture experiments [IH6] , and also extracted indirectly from Coulomb dissociation 
data [TUTU] . Interpretation of Coulomb dissociation data for the capture rate requires careful 
treatment of the parent 15 C and daughter 14 C nuclei in the strong Coulomb field of a heavy 
nucleus besides the nuclear interactions [TTHT3] . Coulomb dissociation provides an alternate 
method to estimate the direct capture reactions involving radioactive isotopes that are often 
difficult to measure experimentally. The 14 C(n,7) 15 C provides an opportunity to compare 
and contrast the capture rates that are obtained from direct capture measurements and 
Coulomb dissociation data [T2l [H]. Developing theoretical methods for radiative capture 
reactions is important for ongoing experimental efforts, and those planned at FRIB |15j . 

We calculate the radiative capture 14 C(n, 7) 15 C cross section at low-energies using halo 
effective field theory (EFT) [T6l ITT] . This reaction has been calculated before in other 
theoretical formulations such as Refs. [21 H21 [TH1 ITS] . Halo EFT has been used to study 
s-wave alpha-alpha resonance [20J and three-body halo nuclei pTJ. Recently it has been 
used to calculate electromagnetic transitions and transition probability strength in one- 
neutron halo n Be [23], radiative neutron capture on 7 Li [221123], and proton- 7 Li interaction 
in coupled-channel extension [25]. In EFT, the cross section is expressed as an expansion in 
the small ratio of low-energy physics scale Q of interest over the high-energy physics scale 
A that involves short distance physics not relevant at low-energy. EFT provides a model- 
independent framework for calculations whose accuracy can be systematically improved as 
long as there is a clear separation between the energy scales, Q <C A. We consider center-of- 
mass (cm.) energies < 2 MeV, corresponding to momenta p < 60 MeV, that is below the 
threshold for the excited states of 14 C nucleus (or neutron). As such in the EFT, the neutron 
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and 14 C core are treated as inert point-like particles. The ground state of 15 C, identified as 
J 71 " = | + , has a neutron separation energy B of only 1.218 MeV that correspond to a binding 
momenta of 7 = y/2j!B pa 46.21 MeV, where \l is the neutron- 14 C reduced mass. In nuclear 
structure calculations the ground state of 15 C can be considered a single neutron halo bound 
to a 14 C core. Then in the single-particle approximation, it is described as a 2 S\/2 state of 
n + 14 C. We use the spectroscopic notation 2S+1 L j with S the spin, L the orbital angular 
momentum and J the total angular momentum. The momenta p, 7 are the soft scale Q. 
The energy threshold for the excited states of 14 C, pion physics, etc., is identified with the 
hard scale A ~ 100 - 200 MeV. 

At low-energy, the capture from lower partial wave initial states should dominate. How- 
ever, neutron capture from initial s-wave state to the ground state through Ml transition 
is suppressed (at one-body current level) due to the orthogonality of the continuum and 
bound state wave functions. The lowest multipole transition to the ground state is through 
El transition from the initial p-wave states 2 P\/2 and 2 P^/2- We note that transition from 
the initial s- and p-wave states to the excited state of 15 C J T = | + is possible. However, 
transitions to the excited state has been found to be a small contribution to the total capture 
rate [21 E3 I2B]- We ignore such contributions in this calculation where we concentrate on 
the dominant effects. 

The paper is organized as follows. In section [IT] we introduce the basic theory and the 
interactions necessary for the 1A C(n, 7) 15 C cross section calculation. The Lagrangian for 
the s- and p-wave interaction of neutron and carbon- 14 is presented. We describe how the 
EFT couplings can be constrained from data. The El capture cross section is calculated in 



section |III| We consider both direct capture and Coulomb dissociation data. EFT couplings 
are constrained to reproduce the available data. From the analysis, we formulate a power 
counting for estimating the sizes of the couplings and the various EFT contributions. In 



section IV we present our conclusions. 



II. FORMALISM 



The construction of the EFT for 14 C(n, 7) 15 C require description of the n+ 14 C bound 
state in the 2 Si/ 2 channel, and the initial state interaction of n+ u C in the 2 Pi/2 and 2 P3/2 
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channels. The interaction in the 2 Si/2 channel is written as 

V 2 



2M 



a + ^ {O) K(^aC)+h.C.], (1) 



where a is an auxiliary field with a spin index a, N a is the neutron field and C is the 
carbon-14 scalar field. M = M n + M c with neutron mass M n = 939.6 MeV and 14 C core 
mass M c = 13044 MeV. Using the equation of motion for the field, it can be integrated out 
of the theory in Eq. 0, and the interaction Lagrangian written entirely in terms of four- 
particle neutron carbon-14 interactions. The non-relativistic s-wave amplitude is calculated 
from the diagrams in Fig. [TJ We get 



where the dressed (f) propagator is 

i 

iDt(p ,p) = A(0)+po _ p 2 /{2M)+t[hm] 2 fo{p ^ p y ( 3 ) 



fo(po,P) = - ftp 



X\ A ~ D f d D ~ l q 1 



2 J J (2tt) d - 1 q 2 - 2pp + pp 2 /M - i0- 



— — ^-(A - V-2/iPo + W 2 /M - z0+), 

Z7T 

with A ~ Q the renormalization scale. We use the power divergence subtraction scheme 
where divergences in space-time dimensions D = 4 and lower are subtracted [27] . In Eq. (j2| , 
we iterate the interaction to all order to describe a s-wave bound state. At low energy 
matching the EFT amplitude Eq. Q to the effective range expansion (ERE) 

■ a , \ Zir i 2-k i 

p cot o — tp p —7 + p[p z + 7 2 J J 2 — ip 

we get 

27rA(0) 1 2 / e x 

M^)P +A=7 -2 P7 ' (5) 

2tt 

where \i = M n M c /(M n + M c ) is the reduced ixiass ; *y ~ 46.21 MeV is the 15 C ground 
state binding momentum and p is the effective range in s-wave. There is no experimental 
constraint on the value of p. A priori it is not clear if the effective range p, which has the 
dimension of length, should scale with the short distance (high-energy) scale p ~ 1/A or 




FIG. 1. Elastic scattering amplitudes AS K ' in s- and p-waves. Double line is the 14 C propagator, 
single line the neutron propagator, dashed line the bare dimer propagator, k = 0, 1, 2 corresponds 
to 2 Si/2, 2 P\/2 an d 2 Pz/2 channels, respectively. 

with the long distance (low-energy) scale p ~ If its the former, p is a next-to-leading 

order (NLO) correction whereas if its the latter, its a leading order (LO) contribution in 
EFT. 

To describe the incoming 2 P\/2 and 2 P-i/2 states we consider a Galilean invariant form 
consisting of the relative neutron and 14 C core velocity vc — v^, and the neutron field N a 
and the scalar carbon-14 field C . In particular we want to project a generic tensor ipf with 
a vector index % = 1,2,3 for the p-wave and a spin index a = 1, 2 for the neutron spin into 
the total angular momentum J = 1/2 piece and J = 3/2 piece. This can be done as 

1 , 



3 ' 



(6) 



where the two pieces are the irreducible forms representing the 2 P\/2 and 2 P%/2 states respec- 
tively. o"j are the Pauli matrices. Thus the p-wave interaction in the EFT can be written 

as 

V 2 



C = y aM 



+ id + 



2M 



X^ + Vsh^lx^P^K 





C + h. cl 



(7) 



where rj = 1, 2 corresponds to the 2 Pi/2 and 2 Ps/2 channels respectively. These particular 
p-wave channels in lx Be were also studied in Refs. [24]. The auxiliary field xf plays in p-wave 



a role similar to the Q field in s-wave earlier in Eq. (jlj). The projectors P°f ,r> in Eq. (|7|) are 



pa 



(8) 



The p-wave elastic scattering amplitude is given by a set of diagrams similar to the s-wave 



amplitude, Fig. [Tj We get 

.;(p7( 2a i),o) = — 



= -[^] 2 ^^(p 2 /(2 /U ),0) = — f- ^ , (9) 



using the p-wave propagator for the x v field 



/ 



"%(po.p)= — — 9lh(v)]2 „, : , (io) 

2^ 



/ l(P », P )= S IC 3 -^A + -A=), 



where ( = ^—2np + fj,p 2 /M — i0 + . 

The EFT couplings in p-wave can be related to observables by comparing the EFT am- 
plitude Eq. (J9]) to the ERE as done for s-wave earlier. For p-wave we get 

■ AT)/ s - 271 " V 2 .27T p 2 

lA{{p) =1 — -T— TTT— -5 — — —7 3 — , (11) 



and 



/i p 3 cot 51 -ip 3 (i _i / a W + r Wp2 / 2 _ 



2 [/^)] 2 2 1 ' 

The ERE parameters , and , rj 2 can in principle be used to determine the 
EFT couplings A*- 1 ), h^ 1 ' and A^ 2 \ in the 2 Pi/2 and the 2 P^/2 channels, respectively. 
However, due to lack of sufficient elastic n+ 14 C scattering data the ERE parameters in 
p-wave are not known. In the EFT it is not clear a priori how the couplings should be 
estimated. In the natural case where all couplings scale with the short-distance scale A, 
initial p-wave interaction would be perturbative. In the presence of shallow bound, virtual 
or resonance states in p-wave, the EFT couplings are fine tuned to scale with powers of 
the long-distance scale Q. Then the p-wave operators in Eq. ([7]) need to be treated non- 
perturbatively [321 [17]. Even in the case where p-wave interaction is perturbative, treating 
it non-perturbatively does not introduce uncontrolled error in the EFT calculation. Thus 
resuming the p-wave interaction with the interactions in Eq. ^ to all order we get a result 
valid in the natural and un-natural case. 

Out of the four unknown p-wave couplings, we can determine two of the couplings from 
the known resonance \ state of 15 C, with a resonance energy E r m 1.885 MeV and width 



T r rj 40 keV in the cm. frame. This resonance state is in the 2 P\/2 channel in the EFT. 
To describe the resonance one needs to treat the p-wave interaction non-perturbatively. 
Analysing the elastic scattering amplitude near the resonance, we get 



This determines the couplings from the resonance parameters. The a\ , ob- 

tained from the \ resonance state when used in the capture cross section Eq. (18) gives 
negligible contribution to 14 C(n, 7) 15 C away from the resonance. Near the resonance it 
produces a sharp peak as we show later in Fig. [3j We determine the scaling of the remaining 
two p-wave EFT couplings by analyzing available 14 C(n, 7) 15 C data in the following. 



III. RESULTS 




(c) 



FIG. 2. El capture. Double dashed line is used to distinguish the final state C dimer field (ft 
from the single dashed dressed dimer field x v representing initial p-wave interaction. Wavy lines 
represent photons, k = 1, 2 corresponds to initial state interaction in the 2 P\/2 and 2 Pz/2 channels, 
respectively. 

The capture reaction lA C(n, 7) 15 C proceeds through the diagrams in Fig. g We only 
concentrate on the El transition. The photon couples to the charge of the 14 C core through 
minimal coupling. This corresponds to gauging the core momentum p — > p + Z c eA, where 
Z c = 6. The contribution from the first diagram Fig. [2] (a) can be projected onto capture 
contribution from initial 2 Pi/2 and 2 Ps/2 channels using the projectors from Eq. (JsJ) - In- 
cluding the contribution from the diagrams (b) and (c) that involve the initial state p-wave 
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interactions from Eq. Q, the amplitude square can be written as 



\M 2p ^\ 2 
9 2P1/2 (P) 



12eh r 



M r 



32M n M c Mp 2 | 2Pi 



/2 



(P)| 



7 — 7 

Air Qti(p 2 + 7 2 



(14) 



P 2 + I 2 -1 /a? + r}V /2 - 

in the 2 -Pi/2 channel. The first term, without the initial state p-wave interaction, in g 2Pl/2 is 
from diagram Fig. [2] (a). In the 2 P3/2 channel we get a similar expression 



\M P ^\ 2 
9 2Pi/2 (p) 



l2eh Q Jz~A 



M, 



16M n M c Mp 2 . 2P . 



+ 



9 

67r/x 



P" + 7" -l/af 3 +riV/2-ip 3 



( 2 )„ 



9 3/2 (P) (5 -3 cos 2 



• 3 s 
7 ip — 7 

47r 67r(p 2 + 7 2 ) 



(15) 



We used cm. kinematics: p the carbon-14 core momentum, k the photon momentum and 
k p = cos 8. There is a contribution from the interference between the two p-wave channels 
that vanish when we average over the angle 9 to calculate the total unpolarized cross section. 
We made the leading order approximation |fc| = ko ~ (p 2 + 7 2 )/(2/i). The wave function 
renormalization factor is related to the residue at the pole of the propagator of the <j) 
particle that represents the 15 C ground state. It is calculated from the dressed cj> propagator 
as 



7-1 



_d_ 
dp Q 



[D<t,(po,p)) 



-1 



n 2 h 2 
27T7 



p =P 2 /(2M)-B 



1 + 



1 



PI 



fn 



(16) 



where B = 7 2 /(2p) rj 1.218 MeV is the ground state binding energy. 
The spin averaged differential cross section in cm. frame is written as 

da 



1 


fc| \M\ 2 


327TS 


P 


2 



(17) 



At LO we can write the Mandelstam variable s ~ (M n + M c ) 2 = M 2 . We write the total 



cross section as 



cr{p) 



1 647ra pj{p 2 + 7 2 ) 



2\g p ^{p)\ 2 + A\g p ^{p)[ 



2 M 2 /j, 2 1 - p 7 
where the electron charge is defined as a = e 2 /(An) = 1/137. 



(18) 



The cross section in Eq. (18) depends on three unknown EFT couplings that can be 
expressed in terms of three ERE parameters: the s-wave effective range p, the 2 P3/2 channel 
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scattering volume and the 2 P^/2 channel "effective range" Written in this form, 

the contributions from Figs. [2] (a), (6) and (c) is model-independent as the ERE parameters 
are not model specific definitions but universal that are in principle directly related to the 
n+ 14 C elastic scattering phase shifts. The total 2 Pi/2 and 2 P3/2 contribution from the tree 
level diagram Fig. [2] (a) without the effective range correction p is around 5/ib. This is 
comparable to the data [El [10] in Fig. [i] but also indicates that effective range p correction 
and/or initial state p-wave interaction is important at LO to explain the data. In the 
natural case ~ 1/A 3 ,r^ ~ A, and initial state p-wave interaction in Fig. [2] (b) and 
(c) is suppressed compared to the diagram (a) by factors of Q 3 /A 3 . Two typical unnatural 
cases in p-wave were considered in Refs. [IS] and [TTJ . In the former ~ 1/Q 3 ,r^ ~ Q 
and the p-wave interaction in all the three diagrams are of the same order. In the latter 
~ 1/(Q 2 A), r| 2 ^ ~ A and the p-wave interaction in diagram (b) and (c) is Q/A suppressed 
compare to diagram (a). We construct a systematic EFT by considering p ~ 1/A and 



,(2) 



1/Q 3 , ~ Q. Then the s-wave effective range p correction is a NLO effect, and 
the 2 Ps/2 interactions are LO. We present only the LO result where the effective range p 
contribution is neglected. 




£cm(MeV) 



FIG. 3. Resonant and non-resonant contribution to El capture cross section cr(£' cm ) in the 2 P\/2 
channel. Solid (blue) curve is the resonant contribution, and dashed (red) curve is the non-resonant 
contribution. 

In the 2 P\/2 channel, the LO cross section is determined by the 15 C ground state binding 
momentum 7, and the | state resonance energy E r and width r r . In Fig. 3J we compare 
the contribution from Fig. [2] (a) to that from Fig. [2] (b), (c). The dashed curve shows the 



non-resonant contribution in the 2 -Pi/ 2 channel and the solid curve shows the | resonant 
contribution (in the same 2 Pi/ 2 channel). As expected the resonant contribution is large 
near the resonance energy E r « 1.885 MeV, and comparatively negligible elsewhere. More 
importantly we notice that the non-resonant contribution is non-negligible throughout the 
energy region. This implies that the interference between the resonant and non-resonant 
contribution in the total cross section is significant as we see later. 




FIG. 4. El capture cross section a{E cm ) with a {2) = -ni/(Q 3 ), rf ] = 2n 2 Q, and Q = 40 MeV. 
Solid (blue) curve uses (ni,n 2 ) = (2, 1.5); dot-dashed (red) curve uses (711,712) = (1.5, 1.2); dashed 
(black) curve uses (711,712) = (0.818, 1.12). Square (maroon) direct capture data from Ref. [6], 
circle (dark blue) Coulomb dissociation data from Ref. [10] . 

In the 2 Ps/2 channel the undetermined ERE parameters are a{\r^ at LO. In Fig. [4] 
we plot the total cross section parametrized by 

«1 2) = -ni/(Q 3 ), r[ 2) = 2n 2 Q for some 
reasonable values of n\ and n 2 of 0(1). We pick Q = 40 MeV. For example, (ni,n 2 ) = (2, 
1.5) and (rii,n 2 ) = (1.5, 1.2) reproduces direct capture data from Ref. [01 . We also show 
Coulomb dissociation data from Ref. [ID] . A %-square fit to the Coulomb dissociation data 
with Q = 40 MeV gives (n\, n 2 ) = (0.818, 1.12). The resonance contribution near E cm w 1.89 
MeV differs from a simple Breit-Wigner form. This is a result of the significant interference 
between the non- resonant and resonant contribution in the 2 Pi/2 channel alluded to earlier 
in discussing Fig. |3j 



Traditionally the cross section a in Eq. 18 is presented in terms of the 5-factor £?, 



(j/y/E cm for use in astrophysical calculation at low-energy |28j. As the capture proceeds 
through p-wave initial states to s-wave final state, the ^-factor is a constant at low-energy [21 
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FIG. 5. El capture S-factor S n = a /^jE cm . We use the same set of parameters (including p = 0) 
and legends as in Fig. |4} 

[2H]- In Fig- [5] we plot the ^-factor S n = a/^E cm using the cross section a from Eq. (18). 
We use the same values of parameters (including p = 0) used in Fig. |4j The three set 
of values for S n at low-energy are consistent within the 30% accuracy expected of the LO 
result. We note that the larger values of S n (solid curve) are close to the values obtained 
in the microscopic calculation in Ref. [19], and the intermediate values of S n (dot-dashed 
curve) are close to the values obtained in the potential model calculation in Ref. [2]. The 
5-factor is a constant at low-energy and expanding it to the lowest order in energy we get 



M c 2 7 (l - fry) 



12 - 4(af } + 2af } ) 7 3 + ([a? J ] 2 + 2[af } ] 2 ) 7 6 + 0(E cm ). (19) 



The contribution from p-wave interaction in the 2 P\/2 channel through a± is negligible 
at low-energy, Fig. |3j The result in Eq. (19) is accurate to NLO at low-energy where 



contributions from p-wave ERE parameters such as r[ , are suppressed. The NLO 
correction to S n at low energy is through the effective range p contribution as seen in 



Eq. (19). 



In Fig. |6j we look at the El reduced transition probability strength [TT], [29] 

dB(El) 9 fiE c 



dE rel 16tt 3 



l -<r(E cm ), (20) 

and compare with available data [TO] . We ignored any recoil and equated E 1 = E rc \ + B. We 
used (rii = 0.818, n<i = 1.12) with Q = 40 MeV. The agreement with data is not surprising 



since the capture cross section in Fig. [4] was extracted using Eq. (20). This assumed negligible 
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nuclear contribution from the Pb target at the forward angles (large impact parameter) in 
Ref. mg. 




£ re i(MeV) 

FIG. 6. B(E1) strength. Solid (blue) curve uses af ) = -0.818/Q 3 , rf ) = 1.12 x 2Q with Q = 40 
MeV. Circle (dark blue) data from Ref. [T0J- 



IV. CONCLUSIONS 

In this work we consider the radiative capture cross section for 14 C(n, 7) 15 C in halo EFT. 
The dominant contribution from El transition between initial p-wave continuum state and 
final s-wave ground state of 15 C is calculated. The EFT is constructed in the single-particle 
approximation taking advantage of the low neutron separation energy in 15 C nuclei. A 
consistent power counting is developed where the leading contribution involve initial state 
p-wave interactions. Both the resonant and non-resonant interaction is considered. 

The EFT result is written in a model-independent form using the ERE parameters. In 
particular, the result depends on the 15 C ground state binding momentum 7, and on the 
scattering parameters a[ , rj 1 and , that encapsulate the interactions in the initial 
2 P\/2 and 2 -P3/2 channels, respectively. The 2 P\/2 parameters are constrained using the 
resonance energy and width of the | resonance state of 15 C. The scattering parameters in 
the 2 i"3/2 channel are estimated from direct capture and Coulomb dissociation data. 

The EFT calculation is shown to be able to describe the energy dependence of the capture 
cross section at the order of the calculation. The EFT couplings constrained from direct cap- 
ture reaction and Coulomb dissociation have values consistent with the EFT power counting. 
The values of the p-wave couplings constrained from the direct capture and Coulomb dissoci- 
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ation data are also consistent with each other within the expected EFT error 0(Q/A) ~ 30% 
on the coupling. The contribution from the resonance in the 2 P\/2 channel differs from a 
simple Breit-Wigner form due to significant interference with the non-resonant contribution 
in this channel. It would be interesting to see if this can be confirmed experimentally with 
more accurate measurements near the resonance energy. Future work should address contri- 
butions from the excited | + state of 15 C to the direct capture reaction 14 C(n, 7) 15 C. Higher 
order contributions from two-body currents should be explored as well. 
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